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Abstract 

We consider two systems of wave equations whose wave-packet solutions have trajectories that 
are altered by the “anomalous velocity” effect of a Berry curvature. The first is the matrix Weyl 
equation describing cyclotron motion of a charged massless fermion. The second is Maxwell equa¬ 
tions for the whispering-gallery modes of light in a cylindrical waveguide. In the case of the massless 
fermion, the anomalous velocity is obscured by the contribution from the magnetic moment. In the 
whispering gallery modes the anomalous velocity causes the circumferential light ray to creep up 
the cylinder at the rate of one wavelength per orbit, and can be identified as a continuous version 
of the Imbert-Federov effect. 

PACS numbers: 
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I. INTRODUCTION 


In many quantum systems the motion of a wave-packet is governed by semiclassical 
equations of the form umi 

0-17 

k = -|^ + e(xxB), (1) 

0-17 

X = + k X O. (2) 

ak 

In the absence of the last term in the second equation these would just be Hamilton’s 
equations for a particle with hamiltonian 'H(x, k) = ^(k) -|- V (x) moving in a magnetic held. 
The additional k x f2 term in ([^ is the anomalous velocity correction to the naive group 
velocity de/d\s.. 

The vector f2 is a function of the kinetic momentum k only, and is a Berry curvature which 
has different origins in different systems. For a Bloch electron in an energy band in a solid 
the curvature accounts for the effects of all other bands. In particle-physics applications the 
curvature arises from the intrinsic angular momentum of the particle. In all cases it affects 
the velocity because different momentum components of a localized wave-packet accumulate 
different geometric phases when both k is changing and the Berry curvature is non-zero [5] . 
These k-dependent geometric phases are just as signihcant in determining the wave-packet 
position as the k-dependent dynamical phases arising from the dispersion equation u = £(k). 

A particularly simple example occurs in in the dynamics of massless relativistic fermions 
and in Weyl semimetals where bands touch at a point. In both systems the wavepackets are 
solutions of a Weyl “half-Dirac” equation and the Berry curvature arises because the spin (or 
pseudo-spin) vector is locked to the direction of the momentum k. The forced precession of 
a spin-F vector is the paradigm in the original Berry-phase paper [6] and the corresponding 
curvature is simply 

n(k) = (3) 

Even this basic example gives rise to much physics — the axial and gauge anomalies [M 
and the chiral magnetic and vortical effect [TUIIT2] . 

The semiclassical analyses that reveal the anomalous velcity in ([^ are quite intricate as 
they have to go beyond the leading WKB ray tracing equation. It is the aim of this paper 
to consider two simple systems in which the predicted anomalous velocity effect can be 
sought directly in the stationary eigenfunctions of underlying wave equation. In both cases 
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the curvature is given by (|^. The hrst (in section 0 is the circular motion of a massless 
charged spin-1/2 particle in a magnetic held. The second (in section III) is the circular 
motion of a spin-1 photon in an optical hbre waveguide. In the hrst case the presence of 
the anomalous velocity is obscured by the coupling of the magnetic held to the particle’s 
magnetic moment. The second unambiguously displays the expected anomalous velocity 
drift. 


II. CYCLOTRON ORBITS 


We start by considering the cyclotron motion of a massless Weyl fermion with positive 
charge e in a magnetic held B = —i?z. The held is derived from a vector potential A = 
B{y, —x)/2 and its downward direction has been chosen so that the particle orbits in an 
anti-clockwise direction about the axis. 

The Weyl Hamiltonian for a right-handed spin-1/2 particle is 

H =-ia ■ {V-ieA), (4) 


where cr = (cxi, ct 2 , as) denotes the Pauli matrices. We are using natural units in which 
h = c = 1, although we will occasionally insert these symbols when it helps to illuminate 
the discussion. 

Acting on functions proportional to we have 

= I + eBLs + + eBa^, (5) 

where = —i{xdy — ydx) = —id^ is the canonical (as opposed to kinetic) angular mo¬ 
mentum and I denotes the 2-by-2 identity matrix. The eigenvalues of the scalar Shrodinger 
operator in parenthesis in ([^ are 

{2n |/| -h / -h 1} + kl, (6) 


and the corresponding eigenfunctions are 

(l«l+l)/2 


^n,l,kAr,<P) = 1 ^ ) 


n\ 


(n+|/|)!" 4 




( 7 ) 


Both n and I are integers and Ln is the associated Laguerre polynomial. When n = 0, 
kz = 0, and I > 0, the wavefunction (/) corresponds to a particle describing a circular 
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cyclotron orbit with the origin as its centre and radius 




( 8 ) 


If we decrease I while staying in the same Landau level (he. by increasing n so as to keep 

I hxed) the classical circular orbit keeps the same radius but its centre moves away from 
the origin and is smeared-out in 6 over the full 27i. When I = 0 the circle passes through 
the origin. For / negative, the energy no longer depends on I and the Landau level keeps n 
hxed while I continues to decrease. The classical orbit still has the original radius, but no 
longer encloses the origin. In particular the case n = kz = 0 and I < 0 corresponds particles 
in the lowest Landau level but with different orbit centres. 

By applying the projection operator P = [E + H)/2E to the Schrodinger eigenfunction 
we hud that the cyclotron-motion eigenfunctions of the Weyl hamiltonian H with n = 0, 
I > 0, and longitudinal momentum kz are 




{Ei^kz + kz)r 
ieB 


exp 


eBP 


( 9 ) 


These states have energy Ei^^z = \/2leB + kl and the orbit radius is still 



At kz = 0, the angular velocity of a wave-packet is 


0 


dl 


fcz=0 


ensuring that = Ri(j) = c = 1. 

There is a special case where I = n = 0 and 



t/’o,o,fcz — e 


ikzZ 



exp 


eB|2p\ 

) 


( 10 ) 


( 11 ) 


( 12 ) 


with E = —kz- This mode only exists as a positive energy mode for kz < 0. It is this 
unbalanced mode, with a density of ei?/27r per unit area in the x, y plane that is the source 
of the chiral-magnetic-effect current 


_ e^B r dkz 
271 Jq 271 




B/i, 


(13) 
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of a gas of zero-temperature Weyl fermions with chemical potential /r mm- 

Consider the I > 0, kz = 0 orbits. Even though these orbits possess no component of 
momentum in the z direction, plugging the time dependence of the classical orbital momen¬ 
tum k into the anomalous-velocity formula (|^ suggests that they should creep down the z 
axis. To compute the predicted creep-rate we observe that a particle of helicity S whose spin 
direction is forced to describe a circle of co-latitude 6* on a sphere with polar co-ordinates 6, 
(j) accumulates Berry phase at the rate |6j 


7Berry = -5'(1 - COS 0)0. 


(14) 


For S = -1-1/2, and using our expression for 0, this becomes 


7Berry = ^(cOS0 - 1)0 = ^(cOS0 - 1)^ 


where 


cos 6 = 


kz 


kz 


( 15 ) 


( 16 ) 


^2leB + kl -MH}' 

In an energy eigenstate this accumulating geometric phase should be indistinguishable from 
the accumulating —Et dynamical phase. In other words yserry should appear as a contribu¬ 
tion to the energy of 


-^Berry = COnst. - - 


kz 


2 ^2leB + kl 



1 kz 

const.--. 

2 21 


This energy adds 


BE 


z = 


Berry 


dkz 


1 1 

2^ 


(17) 


(18) 


to the group velocity over above that expected from the E = c|k| energy-momentum relation. 


The velocity (18) corresponds to drift rate of one-half of a de Broglie wavelength ABrogUe per 
orbit — the de Broglie wavelength being here dehned as the wavelength associated with the 
kinetic momentum k so that E = 27rhc/ABrogiie- 

Unfortunately, except for the special case I = 0, there is no sign of any contribution to 
the energy linear in kz in the exact solution of the eigenvalue problem! Instead we have 


1 


E = \/2leB + kl = V2leB + -k^/V2leB + 


(19) 
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The reason for the absence is that there is another linear-in-fc^ contribution to the energy 
coming from the Weyl particle’s magnetic moment [T^ . 

One might question whether a massless particle can have a magnetic moment as no time 
ever passes on a null-vector world-line. Nonetheless, in the laboratory frame, a particle 
obeying the Weyl equation possesses an energy-dependent effective moment of 

ck 

A^Weyi = positive/negative helicity). 

This moment is precisely what is required for the Larmor precession frequency 


n 


Larmor 


-B 


S 


to coincide with the orbital frequency and so ensure that the spin remains aligned with 
the momentum. The appendix contains a derivation of A^weyi from the Weyl equation and 
explains how the Dirac gyromagnetic ratio g = 2 continues to to be valid even for massless 
particles. For a positive-helicity particle the interaction of this moment with the magnetic 
field provides an energy shift of 

eB ■ k eB cos 9 Ikz 

<3i?magnetic = -^Weyl ' B = — = = +2^' 

that precisely cancels the Berry phase contribution. This cancellation does not apply to all 
eigenstates and does not affect the total chiral magnetic effect, but as explained in [12] it 
does complicate its simple semiclassical derivation in [^. 


III. PHOTONS IN WHISPERING-GALLERY MODES 

We were thwarted in our attempt to observe an anomalous drift velocity in an exact 
solution of the Weyl equation because the accumulating geometric Berry phase was obscured 
by a dynamical phase arising from the particle’s magnetic moment. We therefore seek the 
effect in the motion of a moment-free massless spinning particle. An obvious candidate is the 
spin-1 photon. To make it clear that a photon should obey similar semiclassical equations 
to a massless fermion it helps to rewrite Maxwell’s equations in a form that affords a direct 
comparison with the Weyl equation. The idea for a such a rewriting is apparently due to 
Riemann. The complete set of Maxwell’s equations equations were hrst published in 1865 
and Riemann died in 1866. The rewriting is nonetheless ascribed to him by Heinrich Weber 
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in his edition of Riemann’s lectures published in 1901 [13]. The idea was independently 
discovered by Silberstein [T3] in 1907, and has been recently extensively championed by the 
Bialynicki-Birula’s [U]. The Riemann-Silberstein equations make use the complex-valued 
helds 


where 


= E ± icB 
= EiiZH, 




( 20 ) 


( 21 ) 


is the wave impedance. Provided that Z does not vary with position or time, we can combine 
the two Maxwell “curl” equations as 


idS!^ = 


( 22 ) 


Here Ciocai = 1 /a//^ vary with position. At non-zero frequency, the two Maxwell 
“divergence” equations follow from curl equations and do not need to be separately imposed. 


Once we dehne the spin-1 generator Sj to be the matrix with entries —ieijk, eq. (22) becomes 
a pair of Weyl equations 

ihdt^^ = iciocai (S • (23) 


one for the left-helicity chiral held and one for the right-helicity held, to make the analogy 
with the Dirac-Weyl equation as close as possible, we have inserted an h on both sides of 


(22) so that we can exhibit the equation in terms of the quantum mechanical momentum 


operator p = —iKV. We conclude that as the direction of the wave-momentum vector 
processes, the photon spin is forced to follow it. The photon wavefunction will then acquire 
a geometric Berry phase that is twice as a large as that of the Weyl fermion, and this phase 
must have a similarly-proportioned ehect on the semi-classical particle trajectory. 

The only significant difference between the spin-1/2 Weyl equation and Maxwell equations 
is that the Maxwell “wave-function” obeys a Majorana condition 






(24) 


that indicates that the photon is its own antiparticle. 

To obtain an optical analogue of cyclotron motion, we consider the whispering-gallery 
modes of light in a step-index optical hbre. In whispering-gallery modes a light beam orbits 
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FIG. 1: The separation of the positive helicity circumferential ray (magenta) from the negative 
helicity ray (red) in a cylindrical dielectric waveguide. Wavepackets drift up and down the fiber at 
the rate of one wavelength per turn. 

the fibre circumferentially rather than propagating along its length. We may think of the 
orbit as a closely spaced sequence of total internal reflections off the step discontinuity in the 
refractive index. (See figure[^. A separation of the orbits of left and right circularly polarized 
beams was predicted from the anomalous velocity equations in ng and experimentally 
verified in na. In this section we seek to derive the separation directly from the Maxwell- 
equation eigenmodes. 

Recall that the refractive index n is obtained from the material parameters as n/c = y/ejl 
and, in natural units where c = 1, the local speed of light is 1/n. We take the axis of the 
fibre as the axis and its core to have refractive index n = rii for r < R. The cladding will 
have index n = n 2 for r > R. All fields will have a tacit factor of so our differential 

operators act only on functions of the transverse co-ordinates x, y. We define parameters 7 
by 

7 ^ = ±(epu^ - kl) (25) 

where the sign is chosen so that 7 is real. 

The transverse field components Ej^ and are expressed in terms of the longitudinal 
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components and by 


1 

E_i_ — {fczV _\_Ez ujfi Gz X V ^Hz^ , 

7 

H_l = {kzV±Hz + ueGz X V±Ez} , (26) 

7 

where, for our step fibre (no gradients of e or /r away from the discontinuity) 

- VlHz + {kl - e^iUJ^)Hz = 0 , 

— V\Ez + {kl — eiiiJ^)Ez = 0. (27) 


To be guided we need k^ — efiu^ to be negative in the core, and positive in the cladding. 


The solutions to (27) are then 


Ez = 


Hz = 


HJ;( 7 ir)e*'®, 

r < R, 

i?Kz( 72 r)e*'®, 

r > R, 

CJi( 7 ir)e*'^ 

r < R, 

T)Ki( 72 r)e*'^, 

r > R. 


(28) 


(29) 


Here J/ and K; are the Bessel function and modified Bessel function respectively. 

Below the hbre cutoff frequency the quantity 7 ^ = e/icu^ — k‘1 is always positive. Light 
is no longer completely confined, so the eigen-frequencies have a negative imaginary that 
implies an exponential decay in time. The corresponding eigenfunctions must have outgoing 
waves at inhnity, and so will be of the form 


Ez 


Hz 


HJ/(7ir)e*'®, 

r < R, 

SH{^^(72r)e*'^ 

r > R, 

C'J/(7ir)e*'®, 

r < R, 

T)Hp^(72r)e*^®. 

r > R, 


(30) 


(31) 


Here H^^ is a Hankel function of the first kind. 

Whispering-gallery modes have small kz and large azimuthal quantum number 1. These 
modes are always below the hbre cutoff, and so 7 ^ does not change sign ai r = R. If we 
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consider the waves near the point (x, y, z) = (R, 0 , 0 ) then we have 


^ dE, uydH,\ 

^ 72 V " dr r de J 

_ z (KdE, dH, 

y 72 y r (96' dr 

H = ^ (h ^ 

^ 72 V " dr r de ) 

^ ^ (KdIR_ dE, \ 
y ^2 QQ dr ) 

On substituting the functional form of the solutions for r < R these become 

Ex = 

Ey = ^ ^^Ji(7ir)A + a;/i7iJ'/(7ir)0^ 

Hx = ^{^zliJ\lir)C +^-^:ii{'yir)j^ 

Hy = ^ ^^Ji(7ir)0 - a;e7iJ'z(7ir)A^ , 

and for r > R 


( 32 ) 


(33) 


Ex = ;^(^A;,72Hf)(72r)5-^Hp)(72r)D^ 

Ey = :|(^^HpH72r)5 + n;/i272Hf)(72r)/^^ 

Hx = ;|(^A;.72Hf^(72r)/^ + ^H«(72r)i?j 

Hy=\{ !MH«(72r)D - a;72e2H'«(72r)i?') . (34) 

72 \ ^ / 

The boundary conditions are that the tangential components E^^ Ey and Ely be contin¬ 
uous. The continuity of the normal components of D and B is then ensured by the Maxwell 
“curl” equations. 

If kz = 0, then 7 = noj = and the boundary condition equations break into two 

blocks. The Ez and Hy continuity equations for one block are, respectively, 


Jz(nia;i?)2l 
J' {niUjR)A 


Af\n2uR)B, 


(35) 
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The Hz and Ey continuity equations for the other are 


Ji{niuR)C 

J'{n^uR)C 


B.\^\n2UjR)D, 

.[^}lf\n2uR)D. 


(36) 


There are therefore two families of whispering-gallery modes. The hrst is comprises TM 
modes (transverse when looking along the hbre) that have = 0 (and therefore non-zero 
A, B) with frequencies determined by the eigenvalue equation 


R{niuR) _ nJ^Ii[^\n2UjR) 

V ei J'liriiuR) V 62 }lf\n2uRy 


The second comprises TE modes with Ez = 0 (and therefore non-zero C, D), with 

lET RiniOjR) _ re^B.y\n2UjR) 

V fii Jh (riiuR) V /i2 {n2uR)' 

Equivalently the eigenfrequencies are the zeros of 


(37) 


(38) 


DetTM(a;i?) = \/ ei/piJ'^(nicai?)!!^^ {n20jR) - \/ e2//i2H'^^^ [n2bjR)^i[n\bjK) 
DetTE(a;i?) = \/{niUjR)B.y^{n2UjR) - ^/JI2/~^A'^^\n2UJR)R{nlU}R). (39) 


In general TE and TM modes with the same I are not degenerate. 

Consider hrst the case of rii < n 2 - This is not the situation in a practical optical hbre 
where the core always has a higher refractive index, but there are still low-Q resonant modes 
whose partial conhnement arises because grazing-angle incidence provides strong rehection. 
In this case, as r increases, the Hankel functions are approaching their asymptotic oscillating 
region 


Ey\n2ur) ~ 


in2Lor 


+ ... 


Trn2Uir 

before J/(niu;r) starts to oscillate. For I large and ni/n 2 


(40) 


1 the RHS of equations (37) 


and (38) are slowly varying functions of a; = riiur. They have a very small real part and a 
negative imaginary part. Consequently the eigenfrequencies of both TE and TM modes are 
close to zeros of Ji{x) but we must give x a negative imaginary part to make the determinants 
vanish. The eigenfrequencies are therefore given by 


niUR = Zi^n = - im,n 

where is the n-th zero of Ji(a;) and rji^n is a positive quantity that dihers for TE and TM 
modes. 
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FIG. 2: Plot of the magnitudes of DetTM(wi?) (blue curve) and DetTE(wf?) (brown curve) as a 
function of uR for I = 500. The refractive indices are ni = 1.0, n 2 = 1.1. The green lines are the 
locations of the first few zeros of Ji{niUjR). The locations of the minima are the real parts of R 
times the resonant eigen frequencies. 


Now we introduce a small k^. The pairs of equations are now coupled, but we see that 
the term with kz on the LHS of the Ey and Hy continuity equations multiplies Jz(a;), and 
this quantity is small at resonance. We therefore neglect these terms. If we look at the 
Ex and continuity equations (whose validity was previously enforced by the others) and 
neglect the RHS radiation helds, we hnd that the eigenvalue equations simplify to 


kzniojif'i(niUjR)A ———ifi(niUjR)C = 0 , 

R 

xlojTi^ 

kzniuyi(niOjR)C ^ - - — ifi(niOjR)A = 0 . 

R 


(41) 


As we are very close to zeros of the Bessel function, we can set + niSuR) 


uiRTi(zi^n) and (41) reduce to 


kzA — il 6 uC = 0, 

kzC + ilnl5uA = 0. (42) 


These two equations coincide once we set (A,C) = (l,±mi), which corresponds to right 
and left circularly polarized light in a medium of refractive index ni. They give a frequency 
shift of 

600 = ±-— kz- (43) 

LTli 
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The longitudinal group velocity of a wave packet centered around = 0 s therefore 


du 

dkz nil 


(44) 


Since niOjR ~ I for these modes and niu = k^ is the wavenumber for the light in the core, 
we can write this equation as 


A 


z = 


glass 


nik^R 2'KniR 


(45) 


where Agiass is the wavelength of the light in the hbre core. As the time for one orbit is 
27iRni, we hnd that the rate of drift is one (in glass) wavelength per orbit. This is exactly 
what we expect from the berry phase argument in section [TTj The photons can only make a 
few orbits, however, before they escape or become depolarized due to the TE mode (which, 
looking along the circumferential ray, is a linearly polarized beam with the E held in the 
radial direction) having a shorter lifetime that the TM (which is a linearly polarized beam 
with the E held in the direction). 

In an actual hbre we have ni > n 2 - In this regime the Bessel function J;(nia;r) begins 
to oscillate while the Hankel function li\^\n 2 Ujr) is still almost real and exponentially de¬ 
creasing. In the region of ujR corresponding to total internal rehection, the helds outside 
the glass decay almost to zero as they would for total internal rehection oh a hat interface 
— but eventually the Hankel function begins to oscillate and the helds become outgoing 
radiation. 

We assume that kz is small, so that we can ignore the kz^s in 71 and 72 . The boundary 
condition equations are then 


WJ_ A 

V R eipi ]j Hi ) 

R{niR)A 

(%^— A - i^\l^C] 
\ R eifii V ^ 1 / 

R{rR)C 


ikzl 


H 


( 1 ) 


R ^2^2 ^ 
Hh(72-R) B 


ikzl 1 


H 


( 1 ) 


R e2H2 
Hh(72-R) D 


{^2R)D 

(72 i?) 5 


(46) 


From now on we make the assumption that the impedance does not change from core to 
cladding— I.e. that Hi /This impedence matching condition might be hard 
to engineer, but given that we are seeking a mathematical illustration of the anomalous- 
velocity equation rather than proposing experimental verihcation it is not unreasonable. 
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FIG. 3: Plot of the magnitude of Det(a;ii)/Hp^(n 2 wi?) (blue curve) against uR for an impedance- 
matched fibre with m = 1.1, n 2 = 1 and I = 500. The brown lines are the locations of the first few 
zeros of Ji{niU}R). The three deep minima are resonances in the high-Q total-internal-reflection 
region. The minima become shallower once the angle of incidence of the rays decreases below 
the critical angle, allowing the light to leak out and the eigen-frequencies to acquire a significant 
imaginary part. 

The matching means that the right and left handed Riemann Silberstein fields do not mix. 
It also ensures the degeneracy of the kz = 0 TE and TM modes, both eigenfrequences being 
determined by the same vanishing condition 

Det(a;i?) = i{niuR)}ll^\n 2 UjR) — W^^\n2UjR)3i(niuR) = 0. 

See figure for a plot that locates the zeros of Det(a;i?). 

For non-zero kz the TE and TM pairs of equations are again coupled, but, accepting the 
impedance-matching condition we can decouple the four equations into a different two pairs 
of equations — one for each helicity. We set 

{A,C) = (l,*y^)X+ 

(B,D) = (l,*y^)F+ (47) 

and similarly, with the sign before the i changed, for X_, YL. Using the notation ci ,2 = 
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FIG. 4: Contour plot (blue is lowest value) of the magnitude of Det{LjR)/li^^\n 2 U]R) In the 
complex ujR plane for the same range of parameters as figure\^ The first three zeros lie very close 
to the real axis. They correspond to modes that are strongly confined by total internal reflection. 
The remaining zeros break away from the axis and their larger imaginary parts indicate that the 
light rays can escape because they are below the critical angle. 

the equations for the “+” pair become 

^cjiliiuR/c,) + u^yiiuR/c,)^ X+ = (^^C2^}1 ["\uR/c 2) + u‘^}lf\uR/c2)^ Y+ 

Ji(ojR/ci)X+ = }1 [^\uR/c2 )Y+ (48) 

Dividing the hrst by the second equation and rearranging gives 

2 _ 2 ^^ 2 ( yi{0oR/Cl) Rf\ujR/c2) \ 

^ \R{UJR/Ci) r\^\ujR/c2) J ’ 

When kz = 0 the vanishing of the RHS is the eigenvalue condition. To hnd du/dkz at 
^2 = 0, we therefore need to compute the derivative of the RHS at the points at which it 
vanishes. To do this we make use of the asymptotic formula 
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FIG. 5: Plot of J;'(x)/J;(x) (blue curve) compared with P/x‘^ — 1 (brown curve) for I = 500. The 
brown curve sits on top of the blue curve except for the spikes near the zeros of fi{x). 



FIG. 6: Plot of the real (blue curve) and imaginary parts (green curve) of com¬ 

pared with y^|l — /2/x^l (brown curve) for I = 500. 


Jz(x) 

which is accurate for large x provided that we stay away from places where lfi(x) vanishes 
[see figure . This condition is satisfied at the points of interest [see figure . We may also 
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use the formulae 


Hf*(x) _ fW 
Hp^(x) "" V a :2 


X <1. 


These last two approximations [see hgure imply that 

H'V _ H" _ /H'Y 


0 , 


(50) 


(51) 


which is not quite right, but the derivative on the LHS is 0(1//) in the region of interest 
and can be neglected. We also note that we can evaluate 


3l{x)J Ji{x) \h{x) J 


(52) 


at the unperturbed eigenvalue by exploiting the fact that J[/Jz = H)/H/ at that point. We 
hnd 


I 

R 


C2^)5K 


Cl 

u^R 

Cl 


(IhL.i 

Vn;2i?2 

l\cl-C2^) 


2^2 


(1^ 


Soj 


6u. 


(53) 


Thus 


Su 




(54) 


The opposite circularly-polarized have an opposite frequency shift. We have therefore re¬ 
covered the same drift equation 

doj c 

dkz fill 



that we found for the ni < n 2 hbre. Recall that this drift is at a rate of one wavelength per 
orbit. 


IV. DISCUSSION 

We can compare the mode-expansion rate of drift with that expected from angular mo¬ 
mentum conservation about an axis perpendicular to the reflection interface. For a particle 
of momentum p and helicity S', and treating the orbit as a series of grazing-angle reflections, 
each deflection through an angle Scf) causes a change in the perpendicular spin component of 
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2S'sin(50/2) that must be compensated for by a change in the orbital angular momentum 
of —|p|(5z. Thus |p|i; = —50. Since = c/rigiass we have 


z = -S 


c 

'^glass I P 1.5 


(56) 


In this picture the drift may be understood as a continuous version of the Imbert-Fedorov 
effect [121 Hn] . 

Agreement with the results of the previous section requires us to identify the magnitude 
of the photon momentum |p| with 27rh/Agiass = '^giass^/c, This is the Minkowski expression 
for the momentum of a photon in a medium — as opposed to Abraham’s expression for 
the momentum which places the Ugiass in the denominator. (For a review the Abraham- 
Minkowski momentum controversy see [ 21 ] ■) Minkowski’s momentum is today understood 
to be the pseudo-momentum which is conserved as a result of the homogeneity of the medium 
m- It is the rotational symmetry of the medium about the normal to the core-cladding 
interface that is responsible for the angular momentum conservation, so the appearance of 
the Minkowski momentum is not surprising. 

That a continuous process of reflection and rotation can transport a circularly polarized 
beam of light parallel to itself through an arbitrary distance is related to the fact that a 
continuous sequence of Lorentz boosts and rotations in free space can perform an arbitrary 
“Wigner translation” of a hnite polarized beam |22j . 
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VI. APPENDIX: GORDON DECOMPOSITION, THE WEYL MAGNETIC MO¬ 
MENT, AND g = 2. 

The orginal Gordon decomposition of the Dirac 4-current [22] shows that for any solution 
0 of the massive Dirac equation 


-m)0 = 0, 


(57) 
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the four-current can be expressed as 


pj'y^p) = 

^[Ppv^pj - (y^^)P)) + -5,(^s^"^), 

TT}j 

(58) 

where 



is the Lorentz generator. 

S'"' = 

(59) 


Gordon’s decomposition of the current into a particle number-flux and bound spin contri¬ 
bution clearly requires m ^ 0. There is, however, a version that is valid in both massive and 
massless cases: assume that exp{—iEt} and make use of the Dirac equation 

in the Hamiltonian form 


dti) = —a ■ Vip — im/3ip, 
dtp) = EVp) ■ ex + imp)/3, 

with /d = 7°, a* = 7°7* and so hnd 

j = ^ {p)^Vp) - {Vp)^)p)) + ^(V X S). 

Here 7 = (7^,7^,7^), and 

S = P)^SP) 


(60) 


(61) 


with 




(62) 


so that 

cr 0 
0 cr 

With the particle-number density identihed with p = p)^pj, we can again interpret the first 
term in the decomposition as the current j^ee = epk/E = epv due to particles moving at 
speed V = k/i?. The second term, jbound = (e/-E)V x S is the current due to the gradients in 
the intrinsic magnetic moment density. The magnetic moment itself is found by integrating 
by parts to show that 





^ Jbound d X 


1 

2 


r X 


-VxS 


d^x 


e 

E 


Sd^x. 


(64) 
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For a single massless particle whose spin-1/2 is locked to the direction k of its kinetic 
momentnm this is 

ek 

Mweyl = (65) 

as claimed in section m 

For the both massive and massless case we also have an expression for the momentnm 
density as part of the symmetric Belinfante-Rosenfeld energy-momentum tensor 

( 66 ) 

Using the Dirac equation we evaluate = (i^,P) to find S = E'l|J^^|J, and 

p = 1 (^t(vv.) - X S. (67) 

(If we used the non-symmetric canonical energy-momentum tensor 

- (V>)7''<^). (68) 


we do would not hnd the bound spin-momentum contribution.) 

Again integrating by parts, we recover the spin contribution to the total angular momen¬ 
tum density as 

jrx Qv X d^x = j Sd^x, (69) 

so the division by 2 in the spin contribution to the momentum density is correct. The 
absence of a division by 2 in the formula for the current reflects the g = 2 gyromagnetic 
ratio of the electron. In other words a spin-density gradient is twice as effective at making 
an electric current as it is at contributing to the momentum. 
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